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1. Introduction 

Superfluidity (including superconductivity) has remained at the center of attention of 
low temperature physics since it was discovered at the beginning of last century. By 
now it is well known that both superfluidity in liquids and gases and superconductivity 
in electron systems can all be characterized by the appearance of off-diagonal long-range 
order (ODLRO)^]. In contrast to the diagonal long-range order (DLRO), ODLRO is 
a quantum phenomenon with long-range phase coherence, not describable in classical 
mechanical terms. On the other hand, Yang has pointed out the possibility that ODLRO 
may occur in a solid, coexisting with DLRO P . The possibility for such a supersolid has 
been explored theoretically or numerically in 1970's by a number of authors [21 El IH El El- 
In particular, Leggett ^ has suggested to study superfluidity in a quantum solid in terms 
of non-classical rotational inertia (NCRI). These have stimulated efforts to search for 
superfluidity in solid helium in the following decades [7|. Recent experiments of Chan 
and Kim have shown evidence of NCRI in solid ^He, either confined in porous medium 
[S] or in a bulk 0, with the observed values of the superfluid fraction (SFF) of the order 
0.01. The experimental progress recently drew much attention to revisiting the theory 
of supersolids [Eliniin!. 

In contrast to the efforts concentrated on defects or vacancies, in this paper we shall 
consider a perfect crystal of ^He, where the number, A^, of atoms precisely equals that of 
the sites, and the single-particle density profile has a discrete translation symmetry. We 
shall show that (i) a many-body wave function can be constructed which indeed exhibits 
both DLRO and ODLRO, (ii) in this state short-range correlations due to interatomic 
hard-core repulsion lead to a condensation of zero-point lattice vibrations with long- 
range phase rigidity, (iii) in the presence of an external velocity field, this state shows 
NCRI, and the associated SFF is estimated to be of order 0.01 in a reasonable range 
of parameters. Finally we point out our wave function has features consistent with the 
experiments jHJ E] . 

2. Many-body wave function for a perfect supersolid 

We start with the construction of a wave function that describes a perfect supersolid. 
In a normal crystal individual atoms are localized and oscillate around their equilibrium 
positions, which form a lattice, so that the density profile is periodic, resulting in 
DLRO. For helium atoms, their small mass makes the zero-point oscillations significant. 
Normally the zero-point motion of individual atoms is incoherent, as in the Einstein 
picture of lattice dynamics. This is what happens in normal solid helium (under 
pressure). We propose that at sufficiently low temperature, due to Bose statistics of 
^He atoms and short-range correlations arising from interatomic hard-core repulsion, 
the zero-point motion of individual atoms may become phase-locked and, therefore, 
an ODLRO and phase rigidity will be developed across the whole system. We shall 
not address the question of under precisely what conditions this will happen, but be 
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concentrated on constructing a wave function that demonstrates that the coexistence of 
DLRO and ODLRO is possible in principle. 

We use a localized wave function, 0(r — Rj), to describe the zero-point motion of 
the atom around a lattice point R,. Then the coherent zero-point motion of the atoms 
in the crystal can be described by the following symmetrized product of single particle 
wave functions: 

TV 

5n0(r.-Ri), (1) 

i=l 

with S the symmetrization operator with respect to rj. The symmetrization S in Eq. (1), 
which is absent for an Einstein crystal, incorporates the fact that all atoms are identical 
bosons. The superposition involved in S of permuted products of 0(rj — Rj) takes 
for granted that the zero-point lattice vibrations of individual atoms must be coherent; 
otherwise it would not make sense to superpose the permuted terms. Moreover, as in the 
description of helium superfluid, we should also include a Jastrow factor to incorporate 
short-range pair correlations. Normally the Jastrow factor is taken to be 

I[Jinj) = l[^M-lvirij)}, (2) 

i<j i<j 

with 7 > and the exponents proportional to the interatomic potential v{rij) of the 
Lennard- Jones type, with a hard-core repulsion plus a weak attractive part. Usually 
the product of wave function (P) and (j2I) is used to describe a perfect quantum solid, 
namely 

<iI = sf[<Piv-R,)l[J{r,,). (3) 

i=l i<j 

Provided that the characteristic width, a, of the localized wave packet is much smaller 
than the hard-core size A, the wave function ^ can be approximated by 

. N /n \ 

^o = ^n E</>(r.-R.) Il'^M, (4) 
ViV! \^j=i J 1,^1 

since the Jastrow factor almost annihilates the cross terms in which two atoms are in the 
same localized wave packet at one site. However, when a becomes comparable to A, the 
wave function (0)) incorporates some new features and may exhibit qualitatively different 
behavior than the wave function 0. The above arguments motivate us to propose \l/o 
as our model wave function to describe, at least approximately, a perfect supersolid, 
and we will proceed to show that this wave function allows appreciable Bose-Einstein 
condensation in a density periodic state. 



3. Proof of Bose-Einstein condensation 



From previous experience with ^He superfluid it is known that the Jastrow factor can 
have a Bose-Einstein condensation into the zero-momentum state fSJ El El IHl ITHj . 
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Similarly by expanding the Jastrow factor in Eq. (0)): 



N / ^ N 



^0 = n I ^ E<l^i^^ - I (V^+ ■ ■ ■) , (5) 



i=l 

a non-vanishing zero-mode part {uq ^ 0) would give rise to macroscopic occupation 
in a single particle state, which has a periodic density modulation. (The value of the 
condensate fraction no depends on the microscopic details, and in the superfiuid case it 
was estimated to be of the order 0.01 US].) Now let us proceed to prove that the wave 
function which is of the general form: {u{rij) = 'yv{rij)) 

N N 

^0 ~ n /(r^) n Ar^,) = n /(r.) nexp{-n(r,,)}, (6) 

fc=l i<j k=l i<j 

has a Bose- Einstein condensation on the single-particle state /(r), provided that |/(r)| 
nowhere vanishes and has an upper bound. Assume the same conditions used in 
Ref. |13| IT^ : There exists a positive constant such that the function u{rij) satisfies 
J2l=iu{ris) > —0, for all t, s, ri,r2, ■ ■ ■ satisfying J2i<j<t'^{''^ij) < The proposition is 
the resulting one-particle density matrix possesses ODLRO: 

lini (r|pi|r') = nor(r)/(r'), (7) 

|r— r'|-^oo 

with uq finite and positive, where 

N rJt 

{r\piW) = — l[dr,%{T,v„...)Mr',r2,---) (8) 

VAf i=2 

with Qn the normalization constant of Eq. ((7j) implies that the wave function © 
has a Bose condensation (macroscopic occupation) in the single particle state /(r). 
To prove Eq.©, we notice in the infinite volume limit, 

^0 = lim / drdr' = ^%±i, (9) 
y^oo 1/2 J /*(r)/(r') Qn 

where Ca^+i is defined as 

Civ+i = / J"^"' E dr^nir, r^, ■ ■ ■)v&o(r', r^, ■ ■ ■)■ (10) 
J /*(r)/(r') f4 

Note that the Jastrow functions in Ca^+i are given by 

N 

k=2 2<i<j 

where = r, rf, = r'. The use of the inequality J2l=i u{ris) > —0 allows us to give a 
lower bound for Ca^+i^ 

C,N+l > Qn+1-, (12) 

K 

where A = minM(r) and k = max|/(r)|^. Therefore the condensate fraction ^ has a 
lower bound: 

no > , 13 

Z K 
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where the density n = N/V is fixed in the thermodynamic hmit, and z denotes the hmit 

(i!±iW^. (14) 

Note Qn can be interpreted as the partition function 

(15) 



I n exp I - 5: w ( I r, - I ) + 2 5: In I / (r, 

i=l i<j i 



of a system of classical particles interacting through the two-body potential ksTeffU^r) 
and in an external potential 2k sTeff In |/(r)| with fugacity z. For our wave function 
/(r) = J2j 0(i" " R-i) is periodic, positive and everywhere nonvanishing, so the external 
potential is also finite everywhere. Hence, the thermodynamic limit of this classical 
system exists, and fugacity z is finite if n < Uc, where Hc is the close-packing density. 
In this way we have proved that at T = our many-body wave function (@)) has a finite 
condensate fraction no, leading to ODLRO, Eq.(|7j), a periodic density profile. 

4. Non-Classical Rotational Inertia 

Next we study the response of our perfect supersolid to an external velocity field. 
Suppose that there are N bosonic atoms enclosed in a cylindrical annulus with internal 
radius R and thickness h <^ R. When the cylinder is rotated at a constant angular 
velocity u about its axis, the free energy F{u) measured in the rest frame is of the form 

in 

F{lo) = F(0) + - AF(^), (16) 

where F{0) is the free energy for u = 0, and Iq = NmR^ is the classical rotational 
inertia. The last term is the NCRI, which is related to the SFF a via 

AF{uo) = hooo^a. (17) 

In the rotating frame (with the azimuthal angles ipi ^ ipi + ut), after the gauge 
transformation 

$^exp(^-z^ ' ^ 

the Schrodinger equation reads |1] 



dt 



$. (19) 



Apart from an additive term, Nmuj'^R/2, in energy, $ satisfies the same equation as $o 
in the absence of rotation, however the boundary conditions are changed to 

$(<^. = 0) = exp I -z27r^— J $((^i = 2ti). (20) 

One can therefore conclude that the energy levels are periodic functions of uj. However 
the free-energy is insensitive to the twisted boundary conditions, either when it is in the 
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normal state ^7] or the wave function of all the particles are localized. Only when (at 
least) one extended single-particle state is macroscopically occupied by Nq atoms, the 
SFF is finite and given by 

This idea is similar to that explaining flux quantization in superconductors jIHj. 
For the superfiuid component, we have 



L 



AF(uj) = iVomin / rfx— iV^Pp, (22) 
J 2m 



where sfpe^^ is the order parameter (x|pi|x), and L is the length of the sample. Here as 
a sample calculation, we only consider a one dimensional lattice along the circumference 
x-direction. The one-particle density p has a periodic structure: 

1 ^ 

p{x) = -Y.A{x-R,), (23) 

i=l 

where \Ri — Ri-\\ = d, and the function A, for simplicity, is taken to be a localized 
Gaussian packet: 

,W 1 e.p(4). (24) 

The phase 6{x) is non- uniform, but satisfies 

»(L)=»(0) + ^^. (25) 

Write 6 = 6q + 66, where = CqX with cq = muR/Ti. It is natural to assume that 59 
has a periodic dependence on x with period d. To minimize AF(u;), the phase gradient 
becomes large only in regions with low density. Therefore the SFF is expected to be 
less than the condensate fraction and to increase with decreasing density modulation. 
We apply a simple variational method to prove this. By fixing a global phase we set 

= ,4(?^Uo, (26) 



(27) 





For 59 in the interval [—d/2,d/2\, we take 



59 = Cod 



X f x\^ (X 



'^-d^'A-d) ^''u 



as a trial function, in which the parameters should satisfy the constraint that C5 = 
— 16ci — 4c3. Noting that Noh'^c^/m = UqIquj'^, the SFF a can be expressed as 

; +-.)+3c3(^) +5c,(H) p. (28) 

-d/2 

The right-hand side of Eq. ()28|) is quadratic in ci and C3 and can be easily minimized. 
The resulting SFF is a function of the Lindermann radio a/d. We plot it in Fig. ^ The 



— = A^min / dx 

no ci,c3 



7 



phase 5 as a function of x is shown in Fig. |21for three typical values oi a/d. Our results 
show that the SFF is vanishingly small when a/d < 0.1. When a/d > 0.1, the SFF 
experiences a rapid increase and reaches the condensate fraction no when a/d is around 
0.35. Our result, obtained analytically by a simple variational calculation, is consistent 
with the previous numerical result obtained for three dimensional fee lattice [HI E] • 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 



a/d 

Figure 1. The ratio of SFF a to condensate fraction no as a function of Lindermann's 
ratio a/d. 




x/d 

Figure 2. The phase 9{ in units of cq) as a function of x/d in one period for three 
values of Lindermann's ratio a/d. 



5. Summary 

In summary, we have suggested that in a perfect quantum bosonic solid, under favorable 
conditions, a conspiracy of Bose statistics, significant zero-point fluctuations and short- 
range correlation due to interatomic hard-core repulsion may lead to ODLRO and 
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macroscopic phase coherence, coexisting with crystalhne DLRO and leading to the 
supersohd behavior, e.g. the appearance of NCRI under vessel rotation. We have 
proposed a many-body wave function and proved that indeed both DLRO and ODLRO 
coexist in such a state. The interatomic potential of the Lennard- Jones type, having 
a large hard core and being mostly attractive outside, with favorable parameters is 
arguably credential to choose our state before other non-supersolid states. We note that 
with this mechanism for supersolid the more closely packed the lattice is, the larger 
the SFF a is. Also the presence of a small number of vacancies or impurities (e.g. 
^He) is harmful to condensation of zero-point vibrations, and thus reduces the SFF. 
Observationally, these two features may be exploited as qualitative indication of the 
underlying mechanism suggested here for supersolid behavior. 

The mechanism for supersolid proposed in this note (section |2]) can be understood 
intuitively in the language of path integral. Consider, say, a one dimensional lattice 
of spacing d. \i d — la < A, where a is the amplitude of zero-point vibration, then 
the probability amplitude for two nearby atoms oscillating completely out of phase will 
become negligible due to huge potential energy of the hard-core repulsion. Therefore, it 
is the phase-locked trajectories of oscillating atoms that have lower potential energy and 
dominate the path integral, namely only nearly in-phase zero-point motion of all atoms 
in the crystal are favorable. It should be noticed that these phase-locked trajectories 
are quantum fluctuating pathes describing the ground state in the framework of path 
integral, which should be distinguished from the acoustic phonon excitation. It is these 
phase- locked trajectories that give rise to long-range phase coherence or phase rigidity 
as described by the first factor in Eq.(0]). We call this phenomenon as condensation of 
zero-point lattice vibrations. 

In recent experiments of Chan and Kim jSl IH] that show NCRI of solid ^He, their 
sample is claimed to be ultrahighly pure (with a stated ^He impurity of 0.3 parts of 
per million). Also parameterwise we note ilj^ that in solid ^He the hard core radius is 
a big fraction (> 0.65) of the lattice spacing; and the ratio a/d can be as big as 0.22 
in a variety of situations. As mentioned above, all these features are favorable to our 
suggested mechanism and to get an appreciable SFF. It is arguable that the condensation 
fraction no is of the order 0.01 as estimated in Ref. then in accordance with our 
estimation. Fig. ^ the SFF of our state is appreciable and may reach the order of 
0.01 if a/d is not less than 0.25. (A likely explanation of why Ref. |3] found the SFF 
a < 10~^ is that the estimate used an analogy of the exchange effect for ^He in a normal 
rather than supersolid phase.) Since the SFF depends sensitively on the microscopic 
parameters, to measure the lattice constant and the NCRI at the same time would be 
very desirable. 
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